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Abstract. We use recent lattice QCD outputs to work out the expansion law of the Universe during the 
cosmological quark-hadron transition. To do so, a suitable technique to exploit both pressure and energy 
density data, with the related error bars, is introduced. We also implement suitable techniques to relate the 
T range where lattice outputs are available with lower and higher T's, for which we test suitable expressions. 
We finally compare the cosmological behavior found using lattice data with the one obtainable in the case 
the transition were first order, although not so far from the crossover transition we studied. Differences are 
small to be tested with cosmological data, but the coming of the era of precision cosmology might open a 
channel to inspect the QCD transition through them. 
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1 Introduction 

' In the very early Universe, strongly interacting matter 
. (SIM) was a quark-hadron plasma. Then cosmic expan- 
' sion lowered the temperature T to values 0(100-200 MeV) 
\ and, for a short period, SIM turned into a hadron gas made 
. of TT and mesonic resonances, plus rare nucleons carrying 

■ the cosmic baryon number B. We shall refer to this process 
, as Cosmic Quark-Hadron transition (CQHt). Soon after, 

■ when T shifts below TOtt (pion mass), the only residual 
, SIM shall be made of the baryons needed to carry B. 

CQHt might leave an imprint on today's cosmic ob- 
' servables. This idea was widely explored in the mid-Eigh- 
ties, when the option that CQHt had been a first order 
phase transition was considered likely. In that time lattice 
QCD outputs were already tentatively used to work out 
the cosmic expansion rate and the time dependence of the 
scale factor and several authors also obtained analytical 
determination of the time dependence of the scale factor, 
a(t), during a first order phase transition, occurring close 
to the critical temperature Tc |T] . One of the main ingredi- 
ents to find a right solution of these problems was taking 
into account the presence of the lepton-photon component 
in the CQHt epoch. 

Lattice outputs were considered to this very aim a few 
times and even quite recently (see, e.g., [5]). Often, a tech- 
nique to exploit them was re-invented and a few errors 
occurred recursively, e.g. the neglect of the lepton-photon 
component. 

A topic of this work will then be a general outline of 
the technique and the problems to exploit lattice data to 
cosmological aims. 



A critical issue, concerning CQHt, was however risen 
by Witten [3] in the mid-Eighties, and this made then 
CQHt a hot subject. He showed that, in the case of a first 
order transition, B would tend to remain in the plasma. 
This led him to suggest that quark nuggets could still ex- 
ist in today's Universe, constituting the cosmic Cold Dark 
Matter (CDM) component. The apparent similarity be- 
tween the present densities of the baryon and CDM com- 
ponents could then find a satisfactory explanation. This 
idea has been widely debated and the interest for it never 
completely faded (see, e.g., [4]). 

Other researchers, although keeping to the idea that 
B kept in the plasma, until it existed, suggested that 
the transition would however reach its completion, lead- 
ing however to an inhomogeneous B distribution. Neu- 
tron diffusion would then rapidly smooth out their inho- 
mogeneities. Proton diffusion, instead, requires a coherent 
electron flow and the electron Thomson cross section is 
so large that proton inhomogeneities may last down to 
a temperature ^ 50keV, being therefore able to affect 
Big-Bang Nucelosynthesis (BBN) ^ . This idea triggered 
a large deal of works [6j , of whom we provide a (partial) 
list up to 1990. 

As a matter of fact, a first order phase transition occurs 
in lattice QCD without dynamical quarks [7]. For dynam- 
ical quarks with vanishing mass, the transition is also first 
order. These outputs were given a great importance until 
physical mass quarks could not be treated on the lattice. 
Then, if up and down quarks are light and the strange 
quark is massive, according to phaenomenological results, 
lattice outputs unequivocously yield no confinement phase 
transition, but a smooth crossover. 
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Some doubts can still exist on the nature of the tran- 
sition leading to the breaking of chiral symmetry. Re- 
cent speculation related this transitions to the origin of 
cosmic magnetic fields [8], although a problem exists on 
their scale, which keeps related to the horizon scale at 
t ~ lO^'^sec. 

If CQHt is not a first order phase transition, its rel- 
evance for cosmology is no longer so great. It must also 
be outlined that the recent analysis of CMB, BAO and 
SNIa data [9], within the frame of a ylCDM cosmology, 
has determined a reduced baryon density parameter Ub = 
0.0227 ± 0.0006, in fair agreement with what is predicted 
by homogeneous BEN, in accordance with the observa- 
tional values of ^H/H and '^He/H ratios. 

Some discrepancies however exist with observational 
^He and Li abundances. They are likely to arise from 
marginal errors either in observations or in specific nuclear 
reaction rates. For instance, slight shifts in the rates of the 
^Be{d,p)2'^He or ^_ffe(a, 7) ''iJe reactions ^Oj could put a 
remedy to the '^Li discrepancy. 

The possibility that the observational values of ^H/H 
and '^He/H ratios are an indication of inhomogeneous 
BBN, or derive from other peculiarities as particle decays 
during of after BBN, is still open, but is really second 
choice [TT] . 

On the contrary, in the late Eighties the situation was 
not so strongly constrained; open cosmologies were then 
considered a valid option, while ylCDM was mostly consid- 
ered just as a counter-example in n-body simulations. The 
underlying idea, in the study of inhomogneous BBN was 
that is could relax the constraints on tob, making it com- 
patible with a present overall density parameter flo ^ 0.2, 
so avoiding the need of non-baryonic DM. 

The present situation is however such that, should 
fresh data reopen the option of first order CQHt, we should 
then take care that its features do not produce too strong 
proton inhomogeneities, radically perturbing the theory- 
observation agreement, as far as BBN is concerned. Only 
in this sense, perhaps, BBN remains a constraint to CQHt. 

Let us then also remind that, although lattice QCD in- 
dicates that CQHt was crossover, this is still the outcome 
of a theoretical elaboration. No experimental confirm can 
be soon available and the most direct pattern to explore 
this physics is likely to be just cosmology. With the arrival 
of the era of precision cosmology, however, this might be- 
come a realistic option. 

Besides of trying to provide a general setup on the 
use of lattice outputs in cosmology, this paper therefore 
aims at evaluating the impact that different CQHt pat- 
terns could have on cosmic observables. We shall then also 
consider the possibility that CQHt is a first order phase 
transition, keeping however quite close to recent lattice 
outputs, and will make a comparison between the cosmic 
evolutions in the two cases. 

A further point we need to explore is the connection 
between the very high T regime, where asymptotic free- 
dom is approached, and the T-range where lattice outputs 
are available. We shall show that an analytical expression, 
which can be also seen as a generalization of the histori- 



cal Bag Model [12], is able to fit high-T lattice outputs, 
connecting them with the temperature range where the 
plasma behaves as an ideal gas. The procedure to provide 
this expression and to fit it to data is one of the outputs 
of this work. 

Finally, we shall also need an expression to describe 
SIM in the very low T regime, when it has turned into a 
hadron gas. Although SIM is then just a minor component 
of the Universe, its impact on the cosmic expansion is not 
yet fully negligible. In this paper we shall describe SIM, 
in the form of a hadron gas, by using a state equation in- 
spired to the Hagedorn [13] model, whose parameters will 
be selected in order to connect it with the state equation 
resulting from lattice outputs. 

Back in the Eighties, a large deal of work concerned 
also bubble nucleation, expansion, and coalescence, during 
CQHt (for a review see, e.g., [H]). In the absence of a 
phase transition they are hardly relevant. 

Recent outputs on SIM bulk viscosity during ion col- 
lisions [TS], however, could lead to conjecture some more 
intricate situation. Viscosity is an indication that particle 
reactions proceed too slowly to allow the system to set- 
tle in the equilibrium configurations of an ideal fiuid. A 
tentative data interpretation could then be that reactions 
putting together 3 quarks have a low rate and, perhaps, 
require an intermediate di-quark state. In this case their 
rate would be controlled by the average di-quark concen- 
tration in high-i? quark gluon plasma. 

Should this interpreation of bulk viscosity be correct, 
we expect viscosity to be negligible in the low-B cosmo- 
logical context. 

One might however wonder whether the low rate for 
3q — > B reactions could however mean that B keeps being 
carried by quarks, almost until CQHt is (almost) com- 
plete. This would mean that the residual plasma compo- 
nent would gradually become richer in B, complicating 
the last stages of CQHt, when bulk viscosity could reap- 
pear and even the option of a first order transition might 
reopen. 

In this work we shall not deepen these speculations, 
also because no significant experimental or lattice data can 
be exploited in such analysis. Accordingly, systems will be 
assumed to evolve through thermodynamical equilibrium 
states and no viscosity or heat conduction will be taken 
into account. 

The plan of the paper is as follows: In the next Sec- 
tion we shall briefiy outline some relations obtained from 
thermodynamical considerations in a cosmological con- 
text. Section 3 will just be devoted to show the lattice 
data we shall be using. In Section 4 we shall fit such data 
with a suitable analytical expression, able to reconnect 
them with the asymptotic freedom regime. In Section 5. 
instead, we shall focus on the low-T region and discuss 
how to apply a Hagedorn-like expression to work out en- 
ergy density and pressure. Then, in Section 6, the use of 
interpolating expressions, in the intermediate range, will 
be discussed. Cosmological considerations begin from Sec- 
tion 7, where lattice outputs are actually used to work out 
the cosmic expansion regime. For the sake of comparison. 
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In Section 8 we also deal with a first order transition, ap- 
proaching lattice data. All that will enable us to plot the 
connections between a (the scale factor), t (time) and T 
(temperature); suitable combinations of such parameters 
and their evolution will be shown in Section 9. In Section 
10 we shall gine some examples on the impact that the 
above detailed relations could have on cosmological ob- 
servables, in the era of precision cosmology. Finally, Sec- 
tion 11 is devoted to drawing our conclusions. 



2 Thermodynamics and cosmology 

If we let apart the 3q ^ B reaction, there can be little 
doubts that, around CQHt, relaxation times of particle 
interactions, including neutrinos, are well below the cos- 
mic time. Similarly, the value of the baryon/entropy ratio, 
receiving no substancial contribution after CQHt, allows 
to neglect the chemical potential associated to B. 

Within such context, starting from the thermodynam- 
ical identity dU = -pdV + TdS and setting e = dU/dV, 
a — dS/dV, we have that 



aT = e + p . 



(1) 



Let us then consider the free-energy F — U — TS, such 
that dF = —p dV~S dT. Comparing its second derivatives 
in respect to V and T when taken in different order, one 
has that 

a = dp/dT (2) 
so that the dependence on T of 



-p + T dp/dT 



(3) 



is fully determined, once p{T) is known. If we then define 
<I>^p/T\ E = e/T\ (4) 



the above equation yields 



E{T)=M{T)+T— 



(5) 



If E and <1> depend only on T, this equation can be seen 
also as a differential equation with ^(T) unknown. It is 
then equivalent to the relation T'^E = d{T^(j))/dT and 
can be soon integrated yielding 



1 

T3 



dTT^Eir) , 



(6) 



Tr being a reference temperature. During CQHt, when <P 
and E actually depend just on T, eqs. ([5]) and ^ can be 
used to work out e from p and viceversa. 

Should the transition be first order, however, we have 
two phases. The quark-gluon plasma would have a pres- 
sure Pqg{T), while the hadron gas has a pressure Ph{T) 
and the critical temperature Tc is when pqg[Tc) = Ph{Tc)- 
It must however also be p'qg{Tc) > p'^iTc) (here ' indicates 
differentiation in respect to T) and therefore, according to 
eq. ©, e,g(Te) > eh{T,). 



When Tc is reached, therefore, hadron bubbles must 
nucleate inside the plasma. If this requires no substantial 
supercooling, the Universe then undergoes an expansion 
at constant temperature, while the fraction of each horizon 
occupied by the plasma decreases and the fraction occu- 
pied by hadrons increases. Over large scales, the cosmic 
energy density then reads 



e = egg(l-A) + e,,A + 3<Z',^T4 . 



(7) 



Here A is the fraction of space occupied by the hadron 
gas, while the pressure of the cosmic lepton-photon com- 
ponent, obtainable from 



1.5627 , 



(8) 



depends on the numbers ^bosjer of the boson, fermion rel- 
ativistic spin states in the thermal soup. The above value 
arises from assuming 14.25 effective spin degrees of free- 
dom, i.e. that /j, particles are still fully relativistic. In Ap- 
pendix A we show that, in the T range considered, this is 
acceptable. When aiming at precise quantitative outputs, 
however, the fact that we are in a region where ~ T 
should not be disregarded. 

If supercooling is quite small or infinitesimal, a first 
order phase transition implies quite a small or infinitesimal 
entropy input. In fact, let S = a^a = [a? /T){e + p) be the 
comoving entropy. If T = Tc — 5T (with infinitesimal 5T), 
the average cosmic pressure will be 



P = Pqg{l - \) + PhX + ^i^T^ 



JT 



'.p{T,) - [eqg{T,){l - A) + e„(Te)A]— (9) 



so that, if 5T does not depend on A, 



dp 



ST 



(10) 



Let us consider then the cosmological context, where the 
space-time metric reads 



ds^ =c^di^ -a^{t)d£^ , 



(11) 



di being the infinitesimal comoving length element. Owing 
to the Friedman equation 



d[a^{e+p)] ^ a^dp 



(12) 



it is then 



S^{a^/T)p-S{f/T)= (13) 
= (a^T) (dp/dX) X+a^dp/dT) {t/T)-S{t/T) (14) 

(dots indicate here time differentiation) and the last two 
term at the r.h.s. cancel, owing to eq. ([2]), so that 

S = a^'-^dEApi^ ^ A = [SM - 5.(T.)]^A 

(15) 
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fully vanishes as (5T/Tc ^ . This equation shows that 
only a tiny fraction of the entropy difference between the 
two phases may enrich cosmic entropy. The rest of the 
difference turns into lepton-photon entropy. This is why 
neglecting the latter component is strongly misleading. 

When dealing with a supposed first order transition, 
we shall then consider a negligible supercooling and S to 
be constant. Any other choice would require further as- 
sumptions. 



3 Lattice QCD data 



Lattice QCD deals with the non-perturbative regime for 
the QCD equation of state. Here we are referring to recent 
results obtained from ^£J. They used a {Nt = 6) x 32'^ size 
lattice to calculate the pressure p and the trace anomaly 
e — 3p for 20 values Ti of temperature. The computation is 
performed using physical values for the masses of the up, 
down and strange quarks. No error is given for pi = p{Ti) 
estimates. The trace anomaly, instead, has a significant 
uncertainty which therefore affects their Ci = e{Ti) esti- 
mates. 

These estimates are reported in the Table herebelow 
(T's in units of 100 MeV, the other quantities in units of 
(100 MeV)4). 

From this table the estimates of the coefficients E{T) 
and <?(T) can be deduced; they will be shown in a number 
of Figures hereafter. 







Table 


I 




T 


e 


P 


Ae 


7 


19 


39300 


12600 


353 


5 


56 


13800 


4300 


164.5 


4 


68 


6800 


2070 


95.2 


4 


19 


4350 


1300 


65.3 


3 


67 


2520 


714 


48.02 


3 


24 


1520 


398 


37.92 


2 


82 


859 


196 


10.54 


2 


61 


635 


127 


10.82 


2 


40 


449 


74.1 


6.63 


2 


27 


341 


47.3 


7.08 


2 


13 


238 


26.6 


5.49 


2 


06 


196 


18.6 


6.3 


2 


03 


170 


15.3 


5.66 


1 


96 


113 


10.1 


4.16 


1 


86 


58.7 


5.35 


2.81 


1 


80 


38.7 


3.57 


1.85 


1 


74 


21.4 


2.45 


2.15 


1 


59 


8.88 


.909 


1.2 


1 


45 


6.24 


.315 


1.15 


1 


40 


2.54 


.189 


0.9 



4 The high-T regime 

At high T, we expect E — e/T"* to approach asymptoti- 
cally a value 



^ - 3<^ - i7Ty30)[16 +{7/ 2) N^oiNf la] 



(16) 



Ncoi and Nfia being the number of colors and flavors. The 
above estimates of e and p show that, at T = 719 MeV, it is 
E ~ 0.94 xE {= 15.627, with Ncoi = 3 and Nfia = 3). In 
order to study the CQHt, however, we need an expression 
telling us how E{T) and <?(T) approach E and <P above 
719 MeV. 

The expression we shall propose, however, extends its 
validity down to T ~ 200 MeV, fitting lattice outputs all 
through this T interval, being then suitable to deal also 
with a wide T-range where lattice outputs are available, 
and not only because it provides an easily treatable fit to 
them. 

Let us then assume 



with 



5o = 



T 



{1+52) 



(17) 



(18) 



an expression which, a priori, depends on the successions 

By using eq. ([5]), it is then possible to work out also 
the recursive expressions; 



with 



T6r 



E^3$ + So{S~$) 1 



T6r-{ 



T5o 



(19) 



1 + .5, 



r+l 



(r = 0,l,....) (20) 



These expressions can be truncated at any order n, by 
assuming 5n = 0. 

If we take 6o = and Sq = 4 we recover the MIT Bag- 
model expression with a bag constant B = "PT*. Such 
expression is known not to fit lattice outputs; however, 
the values Sg ^ 1.48 and To — 145 MeV, are the best-fit 
to the highest 6 <P points. The quality of this fit is shown 
in Figure [TJ 

A much better fit is however obtainable if we keep 
So and set 5i — 0. In Figure [1] we show some results 
obtained by freely selecting sq and working then out si, To 
and Ti from the highest 6 <P points, as above. The detailed 
procedure followed in described in Appendix B. 

Clearly, the presence of the 5i term has the result of 
increasing the bending in the middle-T area: the expres- 
sion (I17p sets a T -point where the pressure turns nega- 
tive. (We do not expect such point to be reached. The 
hadron gas should replace the quark-gluon plasma before 
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T 1 1 , 1 ^ 1 , 1 1 1 1 1 1 1 1 r 




T/lOOMeV 

Fig. 1. Lattice outputs compared with the high-T Etting ex- 
pressions. The E{T) points are reported with their 1-a errors. 
The <1>{T) points are reported without errors (not provided 
in the lattice outputs used). Dashed curves are obtained for 

50 = 0. Continuous curves, instead, are obtaned for So and 

51 = 0. They are fits obtained with the procedure described 
in Appendix B, however freely selecting the So parameter. In 
the case considered, the expressions i)J7D . i)J8l) allow to ap- 
proximate up to 10 (over 20) E lattice values (and ~ 14 # 
points). They are no longer suitable when the curvature of 
$(T) changes direction. In the inner frame, the low-T behav- 
ior of <I>(T) is magnihed. The upper black line yields E, i.e. 
E(T) — 3${T) for an ideal gas of relativistic quarks and glu- 
ons. 



then.) Adding a Sq (> 0) correction means that, while T 
decreases (and the quark-quark distance increases), the 
T-point gradually shifts to a higher value. 

One could argue that subtracting from the ideal gas 
coefficient the power law term (To/T)^° is a reasonable 
way to model the binding action due to gluon exchange, 
when the cosmic expansion tries to separe quarks beyond 
their confinement distance. The swiching on of the second 
power law correction 5o = {Ti/Ty^ could then be reminis- 
cent of the expectation that, as the inter-quark distance 
increases furthermore, the exchanged confining gluons, be- 
ing colored, may become source of further gluon emission, 
so causing a further strengthening of their binding action. 

From Figure [T] we also see that reasonable sq values 
are between 1 and 1.1 . Already at sq = 1.1 the E{T) de- 
pendence exhibits a rebounce, which becomes more pro- 
nounced and shifts to higher T's for greater sq values; 
So = 1-1 appears then as an optimal choice to meet nine 
E{T) error bars (for Ti with i = 12, .... 20), with no re- 
bounce above T12. If one takes a smaller sq (e.i?., sq = 1), it 



is possible to meet also the Tn error bar, at the expences of 
shifting below T12 and T13. (Figure [1] apparently indicates 
that the fit with 'P{Ti) values is fair, even down to i ~ 6-7, 
well below 12. In the inner box we however magnify the 
low-T area and show that, although the absolute values 
of (l>(Ti) are really well approached, the expression (flT)) 
misses the shift from positive to negative second deriva- 
tive, essential to avoid the rebounce of E{T) towards low 
T's.) 

As a matter of fact, a reliable analytical expression, 
meeting half lattice points and connecting the range of lat- 
tice data with asymptotic freedom, is to be implemented 
by a numerical fit to lattice data at lower T's. It seems 
clear that the expressions (fTT]) and (fTS]) . with Si = and 
So — 1.1 are suitable to this aim. 

This choice is somehow corroborated by some curious 
numerical regularities. When sq = 1.1 is chosen, the data 
fit returns si ~ 2so = 2.2, To = 76.56 MeV and Ti = 
229.7 MeV; so that the Ti/T^ ratio differs from 3 by - 2 : 
10^. Such regularities hold in a small interval around 1.1 
(typically between 1.08 and 1.12). 

The expressions with so — 0.8 will also be used, when 
aiming to mimic a first-order phase transition. In this ca- 
se, the fitting procedure yields: si = 2.00, To = 26.1 MeV 
and Ti - 447 MeV. 

We also tried to add an extra term to eqs. (fT7|) -p ^ . 
allowing for 61 ^ 0. The fitting procedure becomes then 
even more intricate. A fit with lattice data is obtained with 
the regular values To = 66 MeV, Ti = 2.15 To, T2 = 3.3 To, 
while So = 1, si = 1.5, S2 = 2. This fit, however, is no 
improvement in respect to the one with (5i = : it actually 
meets the E{T) error bars more at their centers, but has 
a rebounce at greater T, making it useless already for 
i — 12; we avoid further graphics complications and we 
refrain from showing it in Figure [1] 



5 The low-T regime 

Expressions of E{T) and <P{T), for T - 100-200 MeV may 
hardly keep any reference to an ideal relativistic gas. As 
a matter of fact, the lightest hadrons are pions and, if 
we tentatively assume for them relativistic ideal gas ex- 
pressions when T ~ m,rj we work out that their mutual 
distance D, in average, is given by 

(C(3) ~ 1.202); this yields D ^ 1.4/to^, a distance compa- 
rable with hadron hard core. Accordingly, hadron proper 
volumes cannot be neglected, let alone hadron-hadron in- 
teractions. Furthermore, as T increases, we can hardly ne- 
glect heavier hadronic resonances. 

There have been several approaches to the thermody- 
namics of the hadronic gas. In any case, however, because 
of their volume and of the temperature, hadron will never 
yield a major cosmic component. Leptons and photons 
will then have a density coefficient i?/^ ~ (7r^/30)15, while 
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hadrons Eg <C (7r^/30)3. Accordingly, for our aims, tak- 
ing the expressions inspired to the Hagedorn model will 
be fully reasonable. Such expressions contain 3 constant, 
that we shall tune to the low-T lattice outputs. 

The pressure of a Hagedorn gas with vanishing chem- 
ical potential reads then 



P{T) 



1 
6^ 



amw(m) / elk z — — — ^ 

Jo exp[E{k,m)/T]^l 

(22) 



with E{k, m) = \Jk'^ + m? and 



I \ (m\' I m 
w{m) = I — exp — 
Vm/ \-^H 



(23) 



Selecting 7 = 5/2 allows a closed expression, using just 
the exponential integral function 



Ei{x) 



dt 



which reads 



p = 'P{T)r with 



<1>{T) = «o {T/THf^ El [mo(T-i - T^')] 



(24) 



and depends on the Hagedorn temperature Th — 200 MeV 
and on the constants ag and mg, related to rh. 

From this expression of (I>{T) and eq. (5) one easily 
works out that 



E{T) = -<?(T) 



ao 



T 



5/2 



exp[mo(T-i -Tg^)] 
Th/T - 1 



(25) 

Figure[2]then shows the <1>{T) and E{T) behaviors at low- 
T, for 5 reasonable choices of the ao and mo constants. 
Although all parameter choices meet the general trend 
and the order of magnitude of lattice outputs, there is a 
specific feature that none of them approaches, the second 
E{T) point, which is however high. 

Accordingly, we shall use such expressions, in associa- 
tion with lattice outputs, just for T < Ti . Almost arbi- 
trarily we chose then ao = 0.28, toq = 600 MeV. The same 
values will also be used when trying to mimic a first order 
transition, up to the temperature where the expression 
of ${T) crosses the high-T expression P7)) . 



6 Fitting lattice outputs 

In order to use lattice outputs in Friedmann equations, 
one needs interpolating them. In principle, mathemati- 
cal and numerical libraries contain excellent interpolating 
routines, namely the splint and spline programs in Numer- 
ical Recipes |18j , making use of a cubic polynomial. 

However, if one uses such recipe to interpolate (P{T) 
(let us also outline that interpolating p{T) yields essen- 
tially identical results, as a counterproof of the outstand- 
ing efficiency of available routines), the critical difficulty 



^ ■ ^ -eooMev / ' ■ ' 

— m(|=500MeV 




1.4 l.( 
T/lOOMeV 



Fig. 2. Lattice outputs compared with the low-T Hagedorn- 
Uke expressions. The E{T) points are reported with their 1- 
a errors. The <f(T) points are reported without errors (not 
provided in the lattice outputs used). Although the expressions 
meet the trend and the order of magnitude of lattice outputs, 
the high level of the secont E{T) point, if real, is clearly missed. 




T/lOOMeV 



Fig. 3. Lattice outputs and their interpolation, as described 
in the text. The cyan solid lines represent the behavior of low— 
T and high-T extrapolating expressions, out of their validity 
range. 
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arises when E{T) is seeken. The passage from <l> to E im- 
plies differentiation, and this blows up small irregularities; 
the resulting E{T) comprises several unmotivated oscilla- 
tions, of clearly numerical origin. 

In particular, there are fair reason to believe that E{T) 
should be always increasing, as also lattice outputs indi- 
cate, towards its asymptotic freedom value. On the con- 
trary, interpolation yields T intervals, between contiguous 
Ti, where E{T) decreases. This effect is particularly evi- 
dent for i ~ 1, 2, 3, ... and at i > 11-12. 

We eliminate this problem with 3 actions: (i) At high T 
we simply use the analytical expressions ([T7| - (fT8|) . start- 
ing from i = 12. (ii) We interpolate the central points of 
E{T), instead of ^{T), and make use of the integral ex- 
pression ^ to obtain (I>{T). This will be done using Ti 
as Tr, and taking the values of E{Ti) and ^(Ti) given by 
Hagedorn-like expressions with our selected parameters, 
(iii) We select the E{Ti) values at j = 2 and 3, at hand 
within a 2-tT error bar, at the minimal distance from the 
central points able to prevent the spline to yield any E{T) 
decrease. Clearly other options are possible, but they ap- 
pear more intricate. 

The overall behaviors of E{T) and ^(T) is then shown 
in Figure [H where the critical low-T region is also mag- 
nified. 




2 4 6 8 10 



T/lOOMeV 

Fig. 4. Lattice outputs approximated by a first order phase 
transition, as described in the text. Notice the vertical line in 
the E{T) curve, where the transition occurs. 



7 Use of lattice outputs in Friedmann 
equations 

The T dependence described in Figure [3] is then used to 
integrate the Friedmann equations. 

It is then convenient to consider first entropy conser- 
vation, yielding that 

[E{T) + ^{T) + A^i^]{aTf = A{3 + •^i^){amf . (26) 

Here a^, Ti are scale factor and temperature at the "ini- 
tial" conditions, where asymptotic freedom holds. In the 
computations here we took Ti = 80 MeV, so to avoid in- 
tereferences with the electroweak transition. From eq. (|26|) 
it is easy to work out the T(a) dependence and the devi- 
ations of the aT product from constant. 

The dynamical Freedman equation then yields 



a 



,E[T{a) 



3<?, 



Tja) 
T 



(27) 



by taking into account that the initial Hubble parameter 
Hi ~ l/2ti (in the radiative expansion regime), it is then 



''da / ^ + ^1^ 
, a \E[T{s)]/Z + ^i^ 



1/2 


■ T, - 







1 (28) 



In this way one obtains a{t) and thence T{t). 

Using this technique we work out the behaviors shown 
in Figures ??-[51 Before discussing them, let us however 
outline how we proceed to build analogous curves holding 
in the case of a phase transition. 



8 Strongly interacting matter in a first-order 
phase transition 

In Figure[4]we show the <P{T) curves we shall use to mimic 
a phase transition: at high-T we shall assume that <?(T) 
is given by an expression (jl7p with 5i = and sq = 0-8; at 
low-T, but also at T > Ti, we shall use the Hagedorn-like 
expression (|24l) with the same parameters as above. The 
two curves are taken above and below the temperature 
Tc where they cross, respectively. Therefore, the overall 
^(T) behavior obtained in this way, although continuous, 
exhibits a (modest) shift of slope at Tc. Accordingly, the 
low- and high-T curves for E{T) do not intersect. In the 
Figure they are connected by a vertical line. 

A cosmological phase transition, occurring close to the 
critical temperature, requires that, after Tc is reached at 
a time th, the Universe stops cooling down. In any volume 
an increasing fraction of space will be then occupied by 
the hadron gas. When the high-T plasma has completely 
vanished, at a time ti, the cosmic temperature restarts 
decreasing. 

During all this process S is conserved, according to 
eq. (15). 

Of course, one might also consider a transition occur- 
ring after a significant supercooling. In order to work out 
the cosmic expansion law, it would be necessary a supple- 
ment of information, yielding e.g., 5T /T as a function of 
A. 

The integration of Friedmann eqs. during a first order 
phase transition occurring without supercooling was first 
performed in [?]. Such integral was then rediscovered a 
few times by various authors. 
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Let pc be the pressure at Tc, including also the lepton- 
photon pressure. Also e and a will include the lepton- 
photon contribution, all through this Section. They will be 
labeled with h (;), to refer to the time when the transition 
begins (ends) and the same labels will be used for time 
and scale factor. Let us also remind that A is the portion 
of space occupied by the low-T phase. 

From the definition of comoving entropy, we have that 



TcS/a 



-Pc 

S/al 



cJhT,{ah/af-p, = eh{l-\)+ei\ (29) 
= {(h +Pc)/Tc. This eq. sets a link be- 



with Uh = S/c^ 
tween a and A 

The dynamical Friedmann eq. reads then 



{d/af = {8TTG/3)[ahT,{ah/a) 



Pc 



(30) 



and, by taking = {a/ ahf j [th/Pc + 1), can be set into 
the form 



du (1 - u2)-i/2 = {67rGpc)^/^dt 
whose integral reads 

u = sm[{6TTGpc)^/^t + A] , 



(31) 



(32) 



A being determined by the conditions at th- 
in the case of the T-dependence of E and shown in 
Figure m we shall use the same equations as in the absence 
of a phase transition for t < th and > ti, and the above 
expressions during the transition. 



9 Scale factor, temperature, time connections 

When the cosmic temperature shifts from some thousands 
MeV to ~ IQOMeV and the CQHt occurs, scale factor and 
time increase by orders of magnitude. 

In order to describe the cosmic evolution, it is then 
worth showing the behavior of precise combinations of a, 
T and t, starting from the initial values a^, Ti and ti. 

Let us then recall first that, in a radiation dominated 
expansion, when the total number of spin states is con- 
stant, the product aT is also constant. Of course, dur- 
ing the CQHt, aT will have robust variations; entropy 
conservation prescribes that, at the end of the transition, 
ajTf = a,Ti{g,/gf f/^ with g,/gf ~ (47.5 + 14.25)/14.25. 
Accordingly, the fact that aT shall increase by a factor 
^ 1.6 can be soon predicted, and is indipendent from the 
detail of the transition. 

Accordingly, in Figures [5H6] we show how aT varies, as 
a function of the scale factor or the temperature, from its 
initial to its final values, however known a priori. 

Let us then notice, in particular, the feature close to 
the transition end, when lattice outputs are used. It is 
due to the anomaly in E(T) at T2 and similar features 
will be present in the next plots. Such features cannot be 
certainly predicted from conservation theorems. 

In Figure[7]we then exhibit the time dependence of T. 
The discrepancies between the two evolutions are evident. 




0.5 
Log(T/100MeV) 



Fig. 5. Dependence of the aT product on tlie temperature. 
The black curve is derived by using lattice data; the cyan curve 
would hold in the case of a phase transition. Notice the verti- 
cal increase of aT while T remains constant during the phase 
transition. 




0.5 1 

Log(a/a.) 



1.5 



Fig. 6. Dependence of the aT product on scale factor. Black 
and cyan colors as in the previous Figure. 
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Log(t/t,) 

Fig. 7. Time dependence of T. Black and cyan colors as in 
previous Figures. Notice again the Bat interval in the cyan 
curve: during a Erst order phase transition T stays constant 
while time elapses. 

during the transition. On the contrary, a slight offset be- 
tween the two curves, after tre transition is completed, is 
just noticeable. 

As a matter of fact, also in this case we must recall 
that, during a radiative expansion when the number of 
spin degrees of freedom does not vary, tT^ should be con- 
stant. This derives from the dynamical Friedmann equa- 
tion, when it is written in the form 

{l/2tf ^ {8Tr^G/9Q)gT^ , (33) 

true if the relation o (x t^^"^ is strictly valid. 

Even though we assume eq. ([55]) to be true before the 
transition, not only during the transition itself eq. 
looses its validity, but also when g would be stabilized at 
a final value, we can only expect that a cx (t+Sty^'^ with a 
suitable St correction. This is what originates the minimal 
offset of the final curves in Figure [T] 

In principle, therefore, the time and scale dependence 
of tT^ are relevant both for the transition period and for 
the final settlement. 

In Figure [5] we show then the explicit dependence of 
tT^ on the scale factor a. Here again the E{T2) anomaly 
causes a clear feature. But the offset at the end of the 
transition, once again, is just barely noticeable. 

This leads us to conclude that, even when we assume 
significantly different E{T) and <P{T) behaviors, the rela- 
tions among T, a and t exhibit discrepacies mostly during 
the transition. We should select completely awkward - and 
unlikely - E(T) and <P{T) to cause a substantial shift in 
T(t) at the transition end. 




^ QQ II I I I I I I I I I I I I I I L 

0.5 1 1.5 

Log(a/ai) 

Fig. 8. Time dependence of tT^. Black and cyan colors as 
in previous Figures. Notice again the feature arising from the 
E{T2) value. 

We then conclude that, once the transition is com- 
pleted, no significant dependence on the relations among 
T, a and t remains: a regime trend is recovered and, if 
no other quantity was altered during the transition, one 
can hardly expect signals from the relative values of such 
parameters. 



10 Observable quantities 

As an example of possible observables, let us consider non- 
relativistic Majorana spinors, mutually annihilating with 
a cross section CTq, such that the annihilation time r 
{n{aaVT))~^ is O(10~^-10~^) sec, occurring then across 
the quark-hadron transition. Here n is the particle number 
density, and vt is their mutual velocity, in the temperature 
conditions considered. 

Their comoving number density — na^ shall obey 
the equation 

dnc/dt + nl (aaVr) /a^ (34) 

and one must know the time dependences of {cTaVT) and 
a'^, in order to integrate such equation. At large t, such 
integral must yield n — n{adg/a)^, adg being the scale 
factor at a time tdg, the decoupling time. This means that 
the t (X proportionality is normalized so that t/tdg = 
{a/odgf. 

We shall however be cautious, before assuming that 
tdg coincides with r, so that n = nitdg). For instance, 
assuming a purely radiative expansion during the decay 
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stages and (aaVx) = const, one obtains that 



n{t) 



n{tdg) 



a-dg 



l + 2{tdg/T)[l~{tdg/ty/^ 

Therefore, at t ^ tdg, it shall be n(t) — n{adg/a)^ with 



(35) 



n = n{tdg)/{l + 2tdg/T) . 



(36) 



Assuming n = n{tdg) is therefore just a zero-th order ap- 
proximation; on the contrary, taking r = tdg, we see that 
the n/n{tdg) ratio is ~ 0.33 . A different T dependence of 
{(JaVT), or a different t dependence of a, would leave most 
of this treatment unchanged, just modifying the factor 2 in 
front of (tdg/r) in eq. ([55)) . Accordingly, the detailed t de- 
pendence of a (or (avT)) may change the residual amount 
of uninteracting Majorana spinors, possibly constituting 
CDM, by a factor 0(2). 

The point is that cosmology is becoming a precision 
science and a variation by a fraction of percent, in the 
abundance of CDM, may soon be measured. Accordingly, 
a percent approximation in the value of ndg, let alone an 
uncertainty by a factor 0(2), might soon become appre- 
ciable. 

If one aims then at connecting cosmological data with 
microphysical data, such as CTq and its T dependence, 
a precise knowledge of a{t) during the CQHt could be 
strictly required. 



11 Conclusions 

This paper was devoted to updating the treatment of 
the cosmological Quark-Hadron transition, taking into ac- 
count recent lattice outputs. 

Such transition had seemed to bear a major relevance, 
in the mid-Eighties, when most reserachers seriously con- 
sidered the option that it could be a first order phase 
transition. When lattice outputs suggested that it was a 
crossover, its interest faded. 

With the coming of the era of precision cosmology, 
however, a detailed knowledge of expansion law and ther- 
mal history, when the time elapsed since the Big-Bang 
was O(10~®-10~^ sec), could become vital to interconnect 
cosmological observables and microphysical parameters. 

In this work we have seen that the most recent lattice 
data are still hard to be used to work out expansion law 
and thermal history. This is due, first of all, to the uncer- 
tainly on the E{T) behavior, deduced from the trace of the 
stress-energy tensor, still given with wide error bars. Fur- 
thermore, even most recent data concern typically 20 tem- 
perature values Ti and, perhaps, they are still not enough; 
we should need more Ti points. 

For instance, at the low-T end, when strongly interact- 
ing matter can be considered a 7r-|-resonances gas, there 
is a clear anomaly in the E{T) trend, however outlined 
by a single point E{T2). Such anomaly, e.g., vanities the 
use of Hagedorn-like expressions in this T-range. In our 
treatment, we took such anomaly into account and saw 



that it causes evident features also in the cosmological ex- 
pansion law. We however legimately wonder whether such 
features, relying on a single estimate among 20, are to be 
considered safe. 

In the attempt to overcome a part of such difficulties, 
however, we have proposed an expression able to fit al- 
most half of the lattice points in the high-T range, and to 
connect them with the very high T regime. In our opinion 
this expression, which a priori depends on 4 parameters, 
but nicely fits more than 20 e and p values, could also 
allow some insight into the physics of the plasma; as a 
matter of fact, the best fit of the above parameters led to 
values unexpectedly regular, so that the expression might 
prove to be something more than a fitting algorithm. 

The main results of this paper are however well sum- 
marized in the Figures and their captions, that can be 
useful to convey the overall message of this work. 
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Appendix A 

In order to evaluate the energy density and pressure 
due to n particles, which are becoming non-relativistic in 
the T-range considered, one must perform the integrals 



dk 



ew{k) 



./O 



dk 



exp[iy(fc)/T] + 1 ' 
k^W-^{k) 



exp[Vl^(fc)/T] + 1 



{Al) 



where W{k) — ^jk"^ + m^, being the mass of the par- 



Til ■■"M 

tides. If one sets then — ET^, = (PT'^, it is easy to 
see that 



E{T) 



dx Z{m^/Tx) 



exp(a;) 



ex.p[xZ{'m^/Tx)] + 1 ' 



3^(r) = ^ / dxx^Z-\mf,/Tx) 



exp(a;) 



explxZ^mfj^/Tx)] + 1 



with Z = a/1 + {m^/TxY . These integrals can be eas- 
ily performed numerically and yield the results shown in 
Figure [9l 

Accordingly, around 100 MeV, energy density and pres- 
sure of the muon component have not yet substantially 
abandoned their ultrarelativistic values. 

Appendix B 

In order to fit the expression p?)) with data, we follow 
these steps: 

We select a value for Sq and use lattice outputs to 
obtain 

i^,^T,[l~^{X,)/^^/'" . (SI) 

We directly use just the six -ipi values for i ~ 15, ...,20. 
According to the expression ()17p . we expect that, in the 
T15-T20 interval, at least, it is 




1 1.5 
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Fig. 9. Temperature dependence of the energy density and 
pressure of the muon thermal component. The vertical dot- 
ted line indicates the mass; the dashed horizontal line are the 
asymptotic ultrarelativistic values. Clearly, forT ~ m^, energy 
density and pressure have already ^ 90% of their asymptotic 
values. 



(with i,j € 15-20) and 



Tr 



rp~Si 
j 



T- 



where, again, i,j,k,r e 15-20. For any choice of i,j,k,r, 
it is then possible to determine a value of si. There are 
however 15 possible combinations, yielding 15 values si, 
that we can average. Otherwise, we can directly treat the 
15 combinations {Tr''^ -T~''')/{T~'' -T'''^) {i,j,k,re 
15-20) as independent data, seeking the best fit si value. 
For reasonable So values we expect the two procedures to 
yield the same Si; we found consistent outputs (within 
1 : 10^) for So between 0.8 and 1.25 . 

Using such Si, we also average among ip values obtain- 
able from eq. (B3), and then eq. (B2) yields a set of To 
values which are also suitably averaged. 

Although using just the 6 top pressure values in this 
procedure, the agreement with many more pressure and 
energy density points eventually follows (as a matter of 
fact, only 4 points would be barely sufficient if we ex- 
clude Tn). 



V'(T) =To + ifT- 



for suitable values of the parameters Tg, si and Lp (= T^*^ ). 
This expression implies that 



(53) 



T \ ^ ^ ^ \ ^ ^ ^ ^ ^ ^ ^ \ ^ ^ r 
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